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Abstract 

We solve a longstanding stability problem for the Kuramoto model of coupled oscilla- 
tors. This system has attracted mathematical attention, in part because of its applications 
in fields ranging from neuroscience to condensed-matter physics, and also because it pro- 
vides a beautiful connection between nonlinear dynamics and statistical mechanics. The 
model consists of a large population of phase oscillators with all-to-all sinusoidal coupling. 
The oscillators' intrinsic frequencies are randomly distributed across the population ac- 
cording to a prescribed probability density, here taken to be unimodal and symmetric 
about its mean. As the coupling between the oscillators is increased, the system sponta- 
neously synchronizes: the oscillators near the center of the frequency distribution lock their 
phases together and run at the same frequency, while those in the tails remain unlocked 
and drift at different frequencies. Although this "partially locked" state has been observed 
in simulations for decades, its stability has never been analyzed mathematically. Part of 
the difficulty is in formulating a reasonable infinite- limit of the model. Here we describe 
such a continuum limit, and prove that the corresponding partially locked state is, in fact, 
neutrally stable, contrary to what one might have expected. The possible implications of 
this result are discussed. 
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1. Introduction 

Collective synchronization occurs throughout the living world, from the rhythmic 
firing of thousands of pacemaker cells in our hearts, to the chorusing of crickets on a 
warm summer evening [Winfree 1967, 1980; Pikovsky et al. 2001; Strogatz 2003]. What 
is remarkable is that these and many other biological populations somehow manage to 
synchronize themselves spontaneously, without any external cue, despite the inevitable 
diversity in the natural frequencies of their constituent oscillators. 

Thirty years ago, Kuramoto introduced an elegant model of such self-synchronizing 
systems [Kuramoto 1975, 1984; for reviews, see Strogatz 2000 and Acebron et al. 2005]. 
Although the model was originally inspired by biology [Winfree 1967], it has since found 
application to many other parts of science and technology. Examples include the mutual 
synchronization of electrochemical oscillators [Kiss et al. 2002], metronomes [Pantaleone 
2002], Josephson junction arrays [Wiesenfeld, Colet and Strogatz 1996], neutrino flavor 
oscillations [Pantaleone 1998], collective atomic recoil lasing [von Cube et al. 2004], audi- 
ences clapping in unison [Neda 2000] , and crowds walking on wobbly footbridges [Strogatz 
et al. 2005]. 

Aside from its scientific applications, the model has also been an object of mathemat- 
ical interest. Its main virtue has always been its tractability. In the limit of an infinite 
number of oscillators, one could "solve the model exactly," in the physicists' sense, as long 
as one was willing to make some plausible assumptions about the stability and conver- 
gence properties of the solutions. Putting these assumptions on a more rigorous footing 
has, however, turned out to be problematic. 

Indeed, Kuramoto himself realized this from the start and was frank about it. For 
instance, in his 1984 book he presents an ingenious formal calculation and then draws 
attention to its limitations. Specifically, he shows that as the coupling between the oscil- 
lators is increased, the zero solution (corresponding to a completely desynchronized state) 
bifurcates supercritically to a nonzero solution (corresponding to a partially synchronized 
state) at a critical value of the coupling strength. He then remarks that the zero solution 
should be stable below threshold and unstable above it, but writes "Surprisingly enough, 
this seemingly obvious fact seems difficult to prove" [Kuramoto 1984, p. 74]. Similarly, he 
points out that the bifurcating solution is expected to be stable above threshold, though 
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"Again, this fact appears to be difficult to prove" [Kuramoto 1984, p. 75]. 

In this paper we settle the second of these issues, namely, the stability of the par- 
tially synchronized state. We find that this state is linearly neutrally stable, rather than 
asymptotically stable. This result may seem puzzling, but there is a precedent for it: the 
same neutral stability was already established fifteen years ago for the zero solution (now 
known as the incoherent state) for coupling strengths below the synchronization threshold 
[Strogatz and MiroUo 1991]. 

The question studied here may be of interest to readers working on stability analyses 
in other parts of nonlinear science, wherever continuity equations arise, such as kinetic 
theory, traffic flow, plasma physics, and fluid mechanics. The problem formulation involves 
a nonlinear partial integro-differential equation, one of whose stationary solutions (the 
partially locked state) contains both a smooth piece and a delta-function piece. To make 
sense of this, we need to work in an appropriate functional-analytic setting, and carry out 
the linear stability analysis in a space of suitable "generalized functions." 

The resulting technical issues are new, at least in this context. They certainly did 
not arise in previous studies of the other stationary states of the Kuramoto model. For 
example, the stability of the incoherent state can be determined with standard methods, at 
both the linear [Strogatz and MiroUo 1991, Crawford 1994] and weakly nonlinear [Bonilla 
et al. 1992, Crawford 1994] levels. The problem is relatively straightforward because the 
incoherent state is described by a smooth (in fact, constant) density of oscillators in phase 
space. The fully locked state is similarly amenable to conventional techniques, as long as 
N is finite. Its stability analysis can be handled with linear algebra [Aeyels and Rogge 
2004, MiroUo and Strogatz 2005] or Lyapunov functions [van Hemmen and Wreszinski 1993, 
Jadbabaie et al. 2004, Chopra and Spong 2006] , since the finite- A?" locked state corresponds 
to a fixed point for an ordinary differential equation. Even the partially locked state is 
susceptible to familiar approaches, if one regularizes the Kuramoto model by adding noise 
terms to it [Sakaguchi 1988]; then the stability of partial locking at onset follows from the 
weakly nonlinear analyses mentioned above [Bonilla et al. 1992, Crawford 1994]. 

But none of these simplifications are available for the problem studied here. Its 
thornier aspects stem from the combination of a continuum limit, the absence of noise, 
the need to work far from the onset of instability, and the singular nature of the partially 
locked state itself. We imagine that a similar mix of ingredients could crop up in stability 
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problems in other parts of nonlinear science, and hence may be of wider interest. 

The goal of this paper is threefold: set up the continuum limit of the Kuramoto 
model in a mathematically precise fashion; describe the fixed states for this model; and 
carry out the linear stability analysis at these states. The third of these is by far the 
most interesting to us, since it has potential to shed light on the still poorly understood 
dynamics of the finite-A/" system [Strogatz 2000, Balmforth and Sassi 2000, Maistrenko 
et al. 2005]. Ultimately, we will achieve a complete understanding of the spectrum of 
the linearized evolution equation for the fixed states of greatest significance, which we call 
special positive fixed states. These are the only candidates for stability; the other stationary 
states turn out to be manifestly unstable. 

The organization follows accordingly. After reviewing the Kuramoto model to es- 
tablish notation (Section 2), we describe its continuum limit (Section 3) and classify its 
corresponding fixed states (Section 4). Included here is the derivation of Kuramoto's 
original self-consistency equation [Kuramoto 1975, 1984], which becomes rigorous in this 
setting. In Section 5 we develop the technical machinery needed to describe the tangent 
space of the model at a fixed state; this is the natural domain for the linear stability 
analysis. We analyze the continuous spectrum of the linearized model in Section 6, and 
derive a characteristic equation whose roots give us the eigenvalues of the linearization in 
Section 7. Finally in Section 8 we prove that the fully locked special positive states are 
linearly stable, but the partially locked special positive states are only neutrally stable, 
since the spectrum contains the entire imaginary axis! The implications of this result for 
the finite-A?" model are far from clear, although this vaguely suggests that one should not 
expect to see any kind of exponential convergence to a stable configuration in the finite-A?" 
system in the range of coupling for which there is only partial locking. 

Before turning to the analytical development, we would like to add a personal note. 
When we began thinking about this stability problem around a decade ago, we found 
ourselves confused by a number of its features. As we had done on an earlier occasion 
[Strogatz et al. 1992], we turned to John David Crawford for advice. J. D. was a brilliant 
mathematical physicist with expertise in bifurcation theory, plasma physics, and nonlin- 
ear science in general. He was also exceptionally generous and a natural teacher. We 
last saw him in the spring of 1998 at a conference on pattern formation at the Institute 
for Mathematics and its Applications in Minneapolis. A few years earlier he had been 
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diagnosed with Burkitt's lymphoma, and when we saw him at the meeting, he was frail 
from chemotherapy but delighted to be able to renew old friendships and to join in the 
scientific discussions. In particular, he became curious about the stability problem that 
is the subject of this paper. The three of us spent a few afternoons working out some 
preliminary calculations. Tragically, J. D. passed away later that summer, at age 44. He 
was very much in our minds as we gathered the fortitude to finish this project, and we're 
sure it would have been completed much sooner had J. D. still been on our team. We are 
honored to dedicate this work to his memory. 

2. The Kuramoto model 

The Kuramoto model is the system 



where N is the number of oscillators, di{t) is the phase of the ith oscillator at time t, ui 
its natural frequency, and K > the coupling strength. The right hand side of (1) defines 
a flow on the A?"- fold torus , which is the natural state space for the system. If the 
frequencies oji have mean oJ, we can go into a moving frame at frequency UJ to transform 
(1) to a system where the frequencies have mean 0. So we can assume U = without loss 
of generality; then fixed points of (1) correspond to phase-locked solutions in the original 
reference frame. We also assume that at least one uji ^ 0; otherwise, (1) is a gradient 
system and is very easy to analyze. 

To characterize the macroscopic state of the system, it is convenient to introduce a 
complex order parameter defined by 



If we think of a state (^i, . . . , ^jv) as an ordered set of N points e*^^ on the unit circle in the 
complex plane, then Re^"^ is just the centroid of this configuration. The radius < -R < 1 
measures the phase coherence of the oscillators and ij) indicates their collective phase. 




(1) 
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Using the order parameter, we can rewrite the governing equations as 

ei=uji + KRsm{i(; -Oi), i=l,...,N. (2) 

For a given set of natural frequencies uii, . . .,(jJn there exists a locking threshold Ki such 
that (1) has fixed points (fully locked states) if and only if K > Kf, furthermore, for 
K > Ki, (1) has a unique stable fixed point up to rotational symmetry, and hence has a 
unique stable fixed point whose order parameter has angle ip = [Aeyels and Rogge 2004, 
MiroUo and Strogatz 2005]. Equation (2) shows that Ki > \uji\ for all z, so Ki will be large 
if just one of the natural frequencies coi is large. So if the natural frequencies Ui are chosen 
randomly with respect to a probability density function (7 on M which has infinite support, 
then as 00 the system (1) will have no fixed points for most selections of ivi. 

Kuramoto's intuition was that one could still predict the asymptotic behavior of the 
system (1) for large N in the absence of fixed points. He guessed that as AT ^ 00 the 
order parameter might still settle down to an almost constant value, despite the incessant 
motion of the unlocked oscillators. Seeking such statistically steady solutions, one can 
assume the order parameter actually is a constant R > and proceed from there. Then 
the oscillators divide into two classes, the locked and drifting oscillators, according to 
whether equation (2) has a fixed point or not; the locked oscillators have natural frequencies 
u)i G [—KR,KR], whereas the drifting oscillators have > KR. We call these kind of 
states partially locked, assuming there are in fact some drifting oscillators (otherwise we say 
the state is fully locked). Kuramoto showed that on average the drifting oscillators make 
no contribution to the order parameter R. Then, by computing the contribution from the 
locked oscillators, he produced a self-consistency equation for R. The A — > 00 limit of 
this equation has a solution i? > if and only if K is larger than some critical coupling 
K(., which Kuramoto computed in terms of the density function g. Numerical simulations 
later confirmed that the size of the order parameter for the system (1) for large N remains 
close to the value of R predicted by Kuramoto's self-consistency equation [Sakaguchi and 
Kuramoto 1986]. 

It's important to understand that Kuramoto's self-consistency equation is only a 
heuristic (albeit deeply insightful) calculation, so unfortunately no precise conclusions 
about the finite-A system can be inferred from it. However, one can introduce an infinite-A 
analogue of Kuramoto's system, which has the advantage that the states analogous to the 
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partially locked configurations described above are fixed states in the infinite- AT model. We 
replace the oscillators 9i and natural frequencies Ui with probability measures pc^, which 
we think of as describing the distribution of the oscillators with natural frequency u) on 
the circle 5"^. Here u ranges over the support of a density function (7, so a state of the 
infinite- Kuramoto model is in effect a family of probability measures parametrized by 
the natural frequencies u. The measures evolve according to an evolution equation mo- 
tivated by the conservation of oscillators; this is a continuity equation, or equivalently, a 
Fokker-Planck equation with no second-order (difi^usion) term. In this setting, Kuramoto's 
heuristic calculation can be made perfectly rigorous. 

3. The Infinite-iV Kuramoto Model 

We now describe the infinite- Kuramoto model. Let O = [—1, 1] or M, and let g{u)) 
be a probability density function on O, which we think of as specifying a distribution of 
natural frequencies. We assume that g{—uj) — g{u))', g is non-increasing on [0, 00) HQ; and 
g is continuous on Q and nonzero on the interior of Q. Two familiar examples are the 
uniform density given by the constant function 1/2 on [—1, 1], and the standard normal 
density function. (For convenience, we extend to be outside Q in the case 0, = [—1,1].) 
As we shall see, these conditions on g are necessary to facilitate many of the calculations 
undertaken in this paper. 

Let Pr{S^) be the space of Borel probability measures on the unit circle. A state for 
the model will be a family G Pr(S'^), parametrized by a; G O. The map OJ ^ puj must 
satisfy at least a mild regularity condition, and to make sense of this we need to put a 
topology on the space Pr{S^). There are various ways of doing this; we choose the one 
best suited to our purposes. 

Consider the Banach space C''{S^) of A;-times continuously differentiable real- valued 
functions on the circle, where A; is a non- negative integer (if A; = then C^{S^) = C{S^), 
the space of continuous functions on S^). The norm on C''{S^) can be defined by 

||0|U = max (|</>(^)| + |0'(^)| + . -.+ 10^(^)1) 

for (j) G C''{S^). We'll be working with the dual spaces C'^{S^)* throughout this paper, 
so it will be helpful to describe their elements as concretely as possible. Any u G C''{S^)* 
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can be represented as follows: 

(0, u)= (j)diJo+ (p'diii + ■■■+ / (p^'^^d/ik 

where //q, • • • , A*A; are signed Borel measures on and is any C'' function on (the 
signed measures fxo, . . . , fXk ^ire not unique). We can express this more succinctly as 

{(f),u) = {l,u) [ (f)dm+ [ (P^'^Ufi (3) 

where m is normalized Lebesgue measure and fi a signed Borel measure with fi{S^) = 
0, which is now uniquely determined by u. The elements of C^{S^)* can be thought 
of as a certain class of generalized functions or distributions on S^, which we call kth- 
order distributions; these are just measures and their first k derivatives, in the sense of 
distributions. We can interpret v — (1, i>)m + {—IYD^^ in this sense. 

In particular, a Borel probability measure ^ on is naturally an element of the dual 
space C^{S^Y ^ with the pairing given by integration: 




This gives an embedding of Pr{S^) in C^{S^)* , and we use the dual norm on C^{S^)* to 
induce a metric on Pr{S^). A distribution z/ G C^{S^Y is a probability measure if and 
only if (l,z/) = 1, and {(j),^) > for any (j) G C^{S^) such that (f> > 0. This shows that 
Pr{S^) is closed in C''{S^)* for aU k > 0. The inclusion map i : C''{S^) C{S^) is a 
compact operator when k > 1, and hence so is its adjoint i* : C{S^)* C''{S^)*. Any 
probability measure has norm 1 when considered as an element of C{S^)*, so Pr{S^) is 
contained in the image of the unit ball under the map i*, and hence Pr{S^) is a compact 
subset of C'iS^)*. 

From now on we insist that k > 1. The compactness of Pr{S^) in C'^{S^)* has some 
desirable consequences. A compact Hausdorff topology cannot be strengthened without 
sacrificing compactness, or weakened without sacrificing the Haiisdorfi' property. This 
implies that the topology on Pr{S^) is the same for all k. Furthermore, the so-called 
weak*-topologies on C'^{S^)* all induce the same topology on Pr{S^). The closure of the 
span of Pr{S^) in C^{S^)* is the subspace of elements u G C''{S^)* that can be represented 
in the form (3) with fi absolutely continuous w.r.t. Lebesgue measure; we call this subspace 
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C'^iS^)*f,s ■ The space C'^{S^)lf^g is separable (it's in fact isomorphic to L^{S^)), and so is 
an ideal choice for a Banach space in which to embed Pr{S^). (The larger Banach spaces 
C'^{S^)* are not separable for all k>0.) 

Now that we have a good topology on Pr{S^), we can officially define the states of the 
infinite- Kuramoto model. We need a regularity condition on the states which will allow 
us to integrate various things; the mildest form of this is the requirement of measurability. 
So we define the states as follows: 

Definition. A state for the inGnite-N Kuramoto model is a measurable map u ^ 
from Q, to Pr{S^). We denote the space of states by S. 

As is the usual practice, we shall identify two states which agree for almost all a; e 1), so 
a state is actually an equivalence class of maps under this relation, but we will usually be 
tacit about this technicality. Since the weak* and dual norm topologies on C''{S^)* induce 
the same topology on Pr{S^), the measurability condition is equivalent to requiring that 
for any C'^ function (j) on S^, the function 



is measurable on Q. The state space S is naturally a closed subset of the Banach space 
L^{fl,C'^{S^)*f^g) of (equivalence classes of) measurable functions from fl to C'^(S'^)*j,g 
that are integrable with respect to the measure g{uj)djj. (See Lang [1993, Chapter VI] for 
background information on integration of functions with values in a Banach space). 

Let > be a constant, which we think of as determining the coupling strength for 
the model. The rest of the ingredients in the Kuramoto model can be defined as follows. 

Definition. Given a state p E S, its order parameter is the complex number 






The vector field associated to p is the function on S 



given by 



v^{e) = LO + KRsm{ilj - 9). 



(5) 



Note that the map 
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is a bounded measurable function of a;, and so the order parameter Re^'^ is well-defined. 
The vector field varies with a;, so it's best to think of this as a family of vector fields v^, 
on parametrized by a; e M, just like the measures P(^. Now we can finally describe the 
equation that drives the infinite- Kuramoto model: 

Definition. The evolution equation for states p E S is 

^^{p^)+D{v^p^)=0. (6) 

Here D is the derivative on S^, interpreted in the sense of distributions: if is a 
function on the circle, then the action of D{vijp^j) on is given by 

{(l),D{v^p^)) ^ - (p'v^dp^. 

The evolution equation is motivated by the conservation of oscillators with a given fre- 
quency u; in fact, (6) is exactly the equation that governs the flow of a measure or distri- 
bution on 5"^ corresponding to the flow determined by the vector held Voj, except for the 
added twist that v^j depends on p through the order parameter, and so (6) is nonlinear. It 
can be shown that there exists a unique solution to (6) for any initial state p, defined for 
all time t. This is not a trivial matter, since the operator D is unbounded on the space 
C''{S^)*, and so the usual existence and uniqueness theorems for ODEs on Banach spaces 
don't apply. We present a brief outline of the proof of this claim, and save the details for 
a future work devoted to the analysis of the nonlinear system defined by (6). 

Let z{t) be a continuous function on [0, T] taking values in the closed unit disc, and 

let 

w^{e) = uj Kljn{z(t)e-'^). 

Replace by tUj^ in (6) to define an uncoupled version of (6). For a given initial state p^, 
we solve this uncoupled equation for each uj to obtain a 1-parameter family t ^ p^^t and 
substitute back in (4) to obtain a new continuous function Z{t) on [0,T]. This defines a 
map ^ : z Z on the space of continuous functions on [0, T] taking values in the unit 
disc. Fixed points of $ correspond to solutions of the original evolution equation (6). For 
T sufficiently small (independent of the initial condition), $ is a contraction map and so 
has a unique fixed point; this proves existence and uniqueness of solutions on [0, T], which 
in turn implies existence and uniqueness for all time t. 
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4. Fixed States 



Our first task is to determine the fixed states for our model. A fixed state p is just a 
solution to the equation 

D{VojPoj) = 

(for almost all a;). The distributions ^ on that satisfy = are constant multiples 
of Lebesgue measure m on S^, which we normalize so that m{S^) = 1. So the fixed states 
satisfy an equation of the form 

VujPuj = C^rn, (7) 

where C^^ is some coefficient function depending on w. If p has order parameter then its 
associated vector field reduces to Vui{0) = oj, so we see that = m for all a;. We call this 
p the incoherent state, as mentioned in the introduction. Note that p does indeed have 
order parameter 0, since 

277 

e'^dO = 





and so the inner integral in (4) is for all u>. Hence the incoherent state is the unique 
fixed state with order parameter 0. 

Now let's try to understand the fixed states which have nonzero order parameter (there 
are a lot of them). If p is such a state, then so is the rotated state p^° given by 

dpt°ie) = dpUe-9o), 

where is any fixed angle. The order parameters for p^° and p are related by the factor 
e*^°, so we can narrow our search to states for which the order parameter has angle ip = 0; 
in other words, we assume the order parameter is some R> 0. 

Definition. The positive fixed states are those fixed states p & S for which the order 
parameter R> 0. 

Plugging in = gives 

v^{e) ^oj - KRsinO. 
Let a; e il. If |a;| > KR, v^{9) ^ for all 9, and equation (7) gives 

dpUd) = — ^^-^de. 

(jO — KR sm 6 
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Because p^; is a probability measure, we must have 



This integral can be evaluated to obtain 

where the ± is the same as the sign of a;; in other words, C^^ is an odd function of u> (this 
observation will be important later). Thus p^^ is completely determined for |a;| > KR: 



2n \u — KRsmO\ 

It's helpful to imagine that for these natural frequencies the oscillators are continuously 
distributed on the circle according to the measure defined above; we call these oscillators 
or frequencies drifting for this reason. By the way, when K is sufficiently large there exist 
positive fixed states which have KR > 1 (we'll construct these in a moment). So when 
Q = [—1, 1], there exist fixed states which have no drifting oscillators; we call these states 
fully locked. 

If |a;| < KR we must have Cj^ = in (7). To see this, observe that away from the 
zeros of v^,, the measure p^, is given by 

dpUO)^ 

u! — KR sm 9 

just as before. But now the denominator changes sign, so cannot be a positive measure 
unless Cu, = 0. Hence p^j must be supported on the roots of v^j. Let 9^; be the root that 
satisfies 

notice that this choice corresponds to the stable fixed point of the one-dimensional flow on 
the circle defined by ^ = v^{9) = u — KRsm9, for |a;| < KR with R regarded as fixed. 
The other root is 9^^ — n — 9^j, and of course it corresponds to the unstable fixed point of 
the fiow on the circle. 

Thus pu, is just a sum of point masses at these two points. Let w{(jj) be the weight of 
the probability measure pu, at the point 9^. Our measurability assumption on the states 
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Pi^ guarantees that is a measurable function on [—KR^ KR] fl f2, taking values in [0, 1]. 
The weight at 9^ is of course 1 — w{ijj). So for |a;| < KR we have 

= {I - w{uj)) 5e^ + w{uj)50*^, (10) 

where we use to denote the unit point mass measure at the point 9. The case we are 
most interested in is when w{uj) = a.e.; in other words, p^ is a unit point mass in the 
right half-plane (so in the general case, w measures the deviation from these preferred 
states). The intuition here is that these special states have the best chance of being stable 
for the full, infinite- A?" system, since their locked oscillators are all located at their stable 
positions, at least with respect to perturbations that don't change the order parameter. 
This observation motivates the following: 

Definition. The special positive Gxed states are those fixed states p e S for which the 
order parameter R> and the weight function w{uj) = almost everywhere. 

If uj is equal to either KR or —KR, then the only probability measure satisfying 
(7) is a unit point mass at ^ or so the values of w at ±KR are irrelevant. Of 
course technically we don't even have to consider this case, since the state p is completely 
determined if we describe p^j for almost all u>. Since the frequencies |a;| < KR have measure 
Pu, concentrated at one or two points, we call these oscillators or frequencies locked. Every 
positive fixed state has some locked oscillators. If all oscillators are locked, we say the 
state p is fully locked; of course, this can only happen when Q = [—1, 1]. Otherwise, we 
call p partially locked. 

The state p satisfies the equation 

r.27r 



which in terms of real and imaginary parts is equivalent to 

R= [ ( [ cos 9 dp^{9)] g{u)dLO and 0^ [ ( [ sin 9 dp^ {9) ] g{u)du. (11) 



JQ, \J0 / JQ, \J0 

We split each of these integrals according to whether |t<;| > KR or < KR. In the first 
case we have 

'\os9dpU9) = r """""ff. , = 
Q Jq u) — KR sm 9 
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since the integrand has a periodic antiderivative on S^. We also have 



• «^ r r ^^^^^^ 

Jo 



u! — KR sin 6 



Jo \uj - KRsine uj J 



de 



KR 
UJ — 2ttC^ 

This last term is an odd function of a;, and so when we integrate it against g{ijj) on the 
set {|<^| > KR} n O we get zero, since g is even. Hence the drifting oscillators make zero 
contribution to the order parameter, in concordance with Kuramoto's original calculation 
[Kuramoto 1975, 1984]. 

The locked oscillators make no contribution to the sine integral since for these fre- 
quencies 

l-2'K 

I sin 6 dpui (0) — sin 9^ = 
Jo 



KR' 

which is an odd function of ou. For the cosine integral we have 



I 



27r 

cosO dpu;{0) = {1 — w{lo)) cos{9ui) + w{uj) cos{9^) 



.KR 

So R and w satisfy the self-consistency equation 



or equivalently 



= j (1 - 2w{KRs)) Vl - g{KRs)ds. 



(12) 



On the other hand if we begin with some R > and a weight function w that satisfies 
(12), then the formulas above define a positive fixed state p with order parameter R, and 
so we have a complete description of the positive fixed states. We summarize these results 
below. 
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Proposition 1. To every Gxed state p e S with order parameter R > we associate a 
measurable weight function w : flf] [—KR, KR] — > [0, 1] that satisfies the self-consistency 
equation (12). Conversely, given R> and a measurable function w : QO [—KR, KR] — > 
[0, 1] which satisfy (12), the state p dehned by equations (9) and (10) is a positive fixed 
state with order parameter R. 

Since most of our arguments treat the locked and drifting frequencies separately, we 
introduce notation for these sets; let 

= nn [-kr, kr] and n'^ = n- nK 

Now let's look in more detail at the special positive fixed states (where w = 0); then (12) 
reduces to the simpler self-consistency equation 

= j \/l - s2 g{KRs)ds. (13) 

We can parametrize all solutions {K,R) to (13) with K, R > as follows. Let t > 0, 
and let p{t) be the state defined by equations (9) and (10) where we substitute t for the 
constant KR. Also define 

f{t) = J a/1 - s2 g{ts)ds = 2 ^ Vl - s2 g{ts)ds, 

and let 

K^f{t)-\ R = tf{t) 

for t e (0, oo). Then KR = t and the state p{t) is a special positive fixed state with order 
parameter R for the model with coupling constant K. (These states are all simply related: 
if we let p = p(l), then we have 

p{t)^ = Pf • 

So the states p{t) are all identical up to a scaling of the frequencies.) 
/ is continuous, positive, and non- increasing on [0, oo). We have 

Vl ds = ^5^, and lim f{t) = 0. 
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/(O) = 2^(0) 



Therefore the image / ((0, oo)) is either ^0, ZI^^^l j or ^0, ^1^21 ^ depending on whether the 
value is taken on at some t > 0. Hence there is a critical coupling constant given by 



7Tg{0) 

such that (13) has solutions if K > Kc, but not if K < Kc- This is essentially Kuramoto's 
derivation of the critical coupling value for his model. 

What happens aX K = K^p. It depends on the behavior of the density function g 
near 0. Let u)q be the largest value of u such that g is constant on [0, a;]. If > (in 
other words, if g is locally constant at 0) then the function / is constant on [0, wq], so 
there is a family of solutions to (13) parametrized by t e (0,a;o] which all have K = Kc 
(Figure 1(b)). However there are no solutions to (13) with K = Kc and > when 
a;o = (Figure l(a),(c)). 

If we rewrite the function R = tf{t) as 



R = 2 



1 - (^) 9{^)duj, 



then we see that i? is a strictly increasing function of t, with image (0, 1) for t > 0. So if 
we plot the parametric curve {K, R) = {f{t)~^, tfit)) i^i the K-R plane, then we obtain a 
curve C in the first quadrant which defines R as an increasing, continuous function of K 
for K > Kc, with perhaps a vertical segment ai K = Kc] R ^ 1 a,s K ^ oo (Figure 1). 

Now suppose we have a positive solution {Kq, Rq) to (12) for a weight function w 
which is not almost everywhere equal to 0. If we set t = KqRq, then (12) shows that 
Kq^ < f(t), so the point (Kq^Rq) will lie on the hyperbola KR = KqRq, in the region 
below the curve C and above the K-axis. Conversely, if {Kq, Rq) is in this region then 
we can construct a positive fixed state with these parameters as follows. Let {K'q, Rq) be 
the point on C that intersects the hyperbola KR = KqRq. Take the special positive fixed 
state with parameters {Kq, Rq) and continuously deform its weight function from w = 
to w — 1/2; the corresponding fixed states' parameters will trace all points below C on 
the hyperbola KR = KqRq. 

To summarize, we see that for each point {K, R) on the curve C there corresponds a 
unique special positive fixed state with parameters K and R; this state always has weight 
function w = 0. And for each point {K, R) in the region between C and the K-axis there 
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exist (actually infinitely many) positive fixed states with those parameters; these states all 
have weight functions that are not a.e. equal to 0. If f2 = [—1, 1], then the points {K,R) 
on or above the hyperbola KR = 1 correspond to fully locked states, whereas points below 
this hyperbola correspond to partially locked states (Figure l(a),(b)). So in this case there 
is a second critical coupling constant 

Ki = f{l)-' 

such that the model has fully locked states if and only if K > Ki. (An equivalent formula 
for the locking threshold Ki was first obtained by Ermentrout [1985].) 

We wish to stress the distinction between Ki and Kc because there seems to be 
occasional confusion about it in the literature. To put it intuitively, suppose that K is 
gradually increased from zero. The system remains completely desynchronized until K 
reaches at which point the first oscillators begin to phase-lock. Thus Kc marks the 
onset of partial locking. With further increases in K, more and more drifting oscillators 
are recruited into the synchronized pack. When K finally reaches Ki, the locking process 
is complete. Now all the oscillators run at the same frequency. 

Hence, partial locking begins at Kc, full locking begins at Ki. Notice that Ki > K^ 
with equality if and only if g is constant on [—1, 1], corresponding to a uniform distribution 
of natural frequencies. If the support of g is M, full locking is never achieved, so it is natural 
to define Ki = oo in this case (Figure 1(c)). 

5. Linearization at Fixed States 

Our next task is to study the linearization of the evolution equation (6) at a fixed state 
p, which we assume is either a positive fixed state or the incoherent state. Our ultimate 
goal is to describe the spectrum of this linearization, which if contained completely in the 
left half plane would establish the asymptotic stability of the fixed state p in the nonlinear 
model. The domain of the linearized model will be the tangent space TpS at p of the state 
space iS, which is a subspace of the Banach space L^{Q, (S^Tabs) ■ recall the relevant 
definitions. 

Definition. Suppose A be a subset of a (real) Banach space E, and p E A. The tangent 
cone TCpA to A at p is the set of x E E for which there exists a function 7 : [0, to) — > A 
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for some to > such that 7(0) = p and 7'(0+) = x; i.e. 



lim 



-X = 0. 



t 



The tangent space to A at p is TpA = TCpA n {-TCpA). 

TCpA is the set of one-sided tangent vectors at p to the set A, and TpA is the set of 
two-sided tangent vectors. TCpA and TpA are always closed subsets of E. If A is convex, 
then TCpA = CpA, where CpA is the convex cone spanned hj A — p. Therefore TCpA is 
closed under addition and multiplication by non-negative scalars, and the tangent space 
TpA = TCpA n {-TCpA) is a closed subspace of E. 

Before tackling the tangent spaces TpS, let's look at the simpler case of just one prob- 
ability measure /j, on S^, and try to understand the tangent space T^Pr{S^) in C^{S^)*^^g 
for some fixed k (it turns out that now we'll need A; > 2 to get everything we want). We can 
explicitly describe these tangent spaces, at least for the two types of measures that occur 
for the fixed states. A (two-sided) tangent vector at /x to Pr{S^) is just the derivative at 
t = of some function 7(t) = in Pr{S^) with ^ = hq and t ranging over some interval 
{—to, to) in M. This derivative, if it exists in C^{S^)*, is the /cth-order distribution rj on 
defined by the rule 



for any C'^ function (J) on S^. For example, let ^{t) = dt be the unit point mass at the point 
t. For k > 2, the derivative of this function at t = is the distribution which assigns to 
any C'^ function (f) its derivative ^'(0); this distribution is just the distributional derivative 
—D5o (we need A; > 2 to insure that 7 is strongly differentiable in the sense described in 
our definition above). 

Pr{S^) is a subset of the hyperplane in C''{S^)*^^g defined by (l,r/) = 1, so any 
r] e T^Pr{S^) must satisfy the linear condition (1,?]) = 0. This is the only constraint on 
the tangent space T^Pr{S^) C C'^{S^)*^^^ if /i is given by dii{6) = a{6)d6 for some smooth 
positive function a on 5"^ (as is the case for the measures p^j for the drifting frequencies 
a;). To see this, let rj be any smooth distribution on 5"^ such that (1, rf) = 0. Then tr] + is 
a probability measure for t is sufficiently small, so 77 e Tp_Pr{S^). The closure of the space 
of smooth distributions u with (I,?/) = in C^{S^)* is the codimension-one subspace of 
C''{S^)*^g defined by (1,?7) = 0, which therefore must be the tangent space T,j_Pr{S^). 
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On the other hand, suppose that ^{t) = jit & Pr{S^) for t e (— to) ^o) has derivative 77 
at t = 0, and jiQ is supported on some compact K <Z S^. If > is a C'^ function which 
vanishes on K'^^ then 

/ > and / (j) d/jQ — 0. 

Taking the derivative at t = shows that (0, 77) = 0. This imphes that the distribution 77 
is when restricted to K^; in other words, r] is also supported on K. Consequently any- 
tangent vector r] e TnPr{S^) satisfies the condition supp(77) C supp(/x). 

The extreme points of the convex set Pr{S^) are the unit point mass measures dp, 
p & S^. As one might expect, the tangent spaces at these points arc fairly small: Ts^Pr{S^) 
is just the one-dimensional space spanned by D6p provided that k > 2 {TspPr{S^) = {0} 
when k — 1, since Ddp ^ (7^(5'^)*^^). To see this, take p = and let be any smooth 
function on with 0(0) = 0'(O) = 0. Suppose ^{t) = G Pr{S^) for t e {-to, to) has 
7(0) = So and 7'(0) = 77. Construct a non-negative smooth function such that 0(0) = 
and 101 < (take 0(^) to be a large multiple of sin^(^/2), for example). Then 

and (0, //o) = 0(0) = 0. (0,77) = by the same argument as above, so (0,77) = as well. 
This, together with the facts that (1,77) = and supp(77) C {0}, imply 77 = cD5o for some 
c. Similarly, if is a linear combination of two distinct unit point masses Sp and Sg, then 
the elements of T,j_Pr{S^) are of the form 

7/ = Co{8p - 5q) + CiD5p + C2D5q 

where the Cj are constants. 

From now on, let's insist that k >2. This insures us a decent supply of tangent vectors 
at the locked states. Now that we understand the tangent spaces for the types of measures 
that occur for the fixed states p e 5, it becomes a relatively straighforward matter to 
describe the tangent spaces TpS we will be working with. The state space S embeds as a 
convex subset of the space L^{^-, C*'^('S'^)*{,g), and we can generalize the arguments above 
to prove the folUowing. 

Proposition 2. The tangent space TpS C L^(ri, C'^(5'-'^)*(,g) at a fixed state p E S consists 
of all 7/ e L^{n, C'iS^)*^^^) such that 7/^ e Tp^Pr{S^) for almost all uj e fl. 
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Proof. In one direction, suppose r] e TpS. Then there exists a map 7 : 1 1— > <S defined on 
some interval {—to, to) with 7(0) = p such that 



hm 



g{uj)duj = 0, 



where the norm is taken in the space C^{S^)* . Let t„ e (0, to) be any sequence converging 
to 0. The (real- valued) functions 



litn)w - Pu 



Vu 



converge to in the space L^(ri, M) with respect to the measure g{u>)du>; hence by passing 
to a subsequence if necessary, we can assume that 

lf{tn)u> Pu 



lim 

n— >oo 







for almost all to & fl. Hence rj^} G Cp^Pr{S'^) = TCp^Pr{S^) for almost all u>. If we choose 
< we get T]^ G —TCp^Pr{S^), so 77^^ G Tp^Pr{S^) for almost all uj. 
Now consider the set V = {rj E L^{Q, C''^(5'"'^)*^g) | i]^ G Tp^Pr{S^) a.e.}; we wish to 
show that V C TpS. Both V and TpS are subspaces of L^{Q, {S^)'^^^^) and TpS is closed, 
so it suffices to prove that TpS contains a set of elements ij e V whose span is dense in 
V. Any r] eV can be expressed as a sum r] = rf + r]^, where r]'',r]^ eV are supported on 
the locked and drifting frequencies respectively, so we can consider these cases separately. 
The most general r] eV supported on the locked frequencies can be expressed as 

for a; G O', where the coefficients Ci are L^-functions of uj w.r.t. the measure g(uj)duj, and 
satisfy the constraints co{lv) = if w{u>) = or 1, C2{u)) = if w{u>) = and ci{u)) = if 
w{u)) = 1. We consider these three terms separately. 

The element rj E V defined by rj^ = co{u!){Sg^ — Sg*^) can be uniformly approximated 
by linear combinations of simpler elements rj eV defined by r]^ = XAi'^)iSeu, ~^o^)j where 
A C 0,'' is measurable {xa denotes the characteristic function of the set A). Because of the 
constraints on the function cq, we can also assume that for some e > 0, e < w{u>) < 1 — e 
for all (jj E A. We must show that this ij E TpS. Define ^{t) by 

l{t),j = Puj + txAiuj){Se^ - S0*J; 
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clearly 7'(0) = rj. Now 7(t)a; is a probability measure as long as |t| < e, so 77 = 7'(0) e TpS. 

The proof for the other coefficients is similar. For the ci-term, it suffices to prove 
that the element ij E V defined by rjuj = xa{^)D50^ is in TpS, where again A C fi' is 
measurable, and this time we assume that for some e > 0, w{u)) < 1 — e for all a; G A. 
Define 7(t) by 

then 7'(0) = —er]. 7(t)a; is a probability measure for all t, and hence 77 = — e~^7'(0) e TpS. 
The argument for the C2-term is exactly the same. Hence any r] E V supported on the 
locked frequencies is in TpS. 

Finally, suppose r] E V is supported on the drifting frequencies, r] is an element 
of the Banach space L^(0'^, C'^(S'^)*^g), taking values in the co dimension-one subspace 
W C C''^('S'^)*^g consisting of distributions orthogonal to the constant function 1 on S^. 
Therefore rj can be uniformly approximated by linear combinations of elements in V of the 
form 

where A C O'^ is measurable and ^ is a smooth distribution in W has no dependence 
on to). We can also assume that A has positive distance from the boundary frequencies 
±KR. It suffices to prove that this rj e TpS. Define ^{t) by 

clearly 7'(0) = 77. Since ^ is a smooth distribution, ^ is a signed measure on given by 
d^{9) = a{6)d6 where a is a smooth function on with integral 0. Now recall that for the 
drifting frequencies, is the measure given by (9); the coefficient function has minimum 
value 

1 ya;2 - {KRf 
2^ |a;| + KR ' 

which is uniformly bounded away from for a; e ^ since we assumed A has positive 
distance from ±KR. The function a is bounded on S^, so if \t\ is sufficiently small, 7(^)0- 
is a probability measure for all a;, and hence r] = 7'(0) e TpS. 

QED 

Strictly speaking TpS depends on k, even though we supress this dependence in the 
notation. As a closed subspace of L^(f2, C'^(S'-'^)*5g), TpS is a Banach space in its own 
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right, and is the natural domain for the hnearization of the evolution equation (6), which 
we turn to next. If we replace p by p + erj in (5), we see that the first-order perturbation 
of the vector field v^, corresponding to a tangent vector r] G TpS is equal to 

/oo 
(sin(r - e),ri^{T))g{uj)duj = K {Sri cos 9 - Cr] sin 6) 
-oo 

(the pairing inside the integral is with respect to the dummy variable r). The coefficients 
Crj and Srj are given by 

Cij — / {cos T,r]^{T))g{Lo)du) and 
Jq 

Sr]= {siiiT, r]u,{T))g{u;)du. 
Jq 

Crj and St] are respectively the perturbations of the real and imaginary parts of the order 
parameter i?e*^ (C for cosine, S for sine). To linearize (6), replace p,^ with p^^ + er]^ and 
gather all the linear terms in e to obtain the equation 

d 

-^iVui) + D(v^r]^ + K {Srj cos 9 - Crj sin 9) p^) = 0. 
This leads to the following definition. 

Definition. The linearized evolution equation at a state p e S is ^{rf) — Lij, where L is 
the linear operator defined by 

{Lr])^ = -D [v^T]^ + K {Srj cos 9 - Crj sin 9) p^) (14) 

on the space TpS. 

The right hand side of (14) is an integrable family of distributions in C''^"'"^(S'^)*^g, but 
not necessarily in C'^(5'^)*^g, since D may map fcth-order distributions to {k + l)st-order 
distributions. Hence the operator L is in general unbounded (the exception, as we shall 
see, is when p is fully locked). However L is a closed, densely-defined operator on the 
Banach space TpS, which is the next best thing to being a bounded operator. 

6. The Spectrum of L for Special Positive States 
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The main goal of this paper is to describe completely the spectrum cr(L) for the 
special positive states. (For all other positive states, the spectrum a{L) contains positive 
numbers, and so these states are linearly unstable, and hence of less interest to us. We'll 
comment more on this at the end of this section.) Since the spectrum may contain complex 
numbers, we need to study the operator L on the complexified tangent space TpS ® C; so 
from here on we will allow the distributions r]^^ to be complex- valued. We'll denote the 
complex Banach space TpS <S> C simply as E. Following tradition we partition cr{L) into 
three parts: its point, continuous and residual spectrum. The point spectrum crp{L) is just 
the set of eigenvalues of L. The continuous spectrum ac{L) is the set of A e C such that 
ker(A/ — L) = {0} and the image Im(A/ — L) is dense in E, but the inverse {XI — L)~^, 
defined on the dense subspace Im(A/ — L), is unbounded. Since L is a closed operator, the 
densely-defined operator (A7 — L)~^ is bounded if and only if it is defined on all of E, or 
equivalently, if and only if the image of XI — L is E; this is a consequence of the closed 
graph theorem (see Kato [1995, p. 166]). So we can also describe ac{L) as the set of A 
such that XI — L is one-to-one, has dense image, but is not surjective. The remainder of 
the spectrum is the residual spectrum (7j.(L), which is therefore the set of A e C such that 
ker(A/ — L) = {0} and Im(A/ — L) is contained in a proper closed subspace of E. In this 
section we begin the analysis of the spectrum of L. 

Henceforth, unless explicitly noted, we assume that p is a special positive state. To 
understand the spectrum of L, it helps to express L = M + B, with the operators M and 
B defined by 

[Bt])^ = -KD{{Sr] cos - Crj sin 9) p^) . 
Notice that M is completely uncoupled, in the sense that {Mrj)^ depends only on r/i^; in 
other words, if we define M^r],^ = —D{v^r]^), then {Mrf}^ = M^r]^. The operator B is 
bounded on E provided that k > 2, which as mentioned earlier is needed to insure that 
the map u ^ Dp^ from Q to C^{S^)* is measurable. B has a codimension-two kernel 
determined by the equations Ct] = Sr] = 0, so the rank of B is only 2; in other words, the 
operators L and M are in a sense very close. The coupling of the oscillators is entirely 
expressed through the operator B, so one can think of the operator M as describing the 
linearized Kuramoto model with the coupling artificially suppressed. 

We can split the tangent space E asa direct sum E = E^®E'^, where E^ and are the 
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subspaces of tangent vectors supported on the locked and drifting frequencies respectively. 
Notice that M preserves both these subspaces. As we saw in the proof of Proposition 2, 
E'' is isomorphic to the space of complex- valued functions L^(n^,C) (w.r.t. the measure 
g{uj)duj on O^). E'^ is the space L^(0'^, W), where W C C''(S'^, C)*^^ is the subspace of 
complex- valued distributions orthogonal to the constant function 1 on 5"^. We can also 
split the tangent space another way, into even and odd subspaces, as follows. The space 
E has an involution defined by the rule 

We call r] even if 77 = ry and odd if rj = —rj. This notion of even and odd behaves much 
like the usual one for functions of one variable: any rj & E can be expressed as a sum of 
an even and odd element in a unique way; if r] is even then the element Dr] is odd, and 
vice versa; and if (j){d) is a smooth function on 5^, then the usual rules for 0?7 apply: if 
is even (in the traditional sense) and 77 is even, then 0?7 is even, etc. Notice that if r] E E 
is even then Srj = 0, since the integrand is an odd function of u>: 

{sine, rj-UO)) = {sine,rju,i-0)) = {sin{-e) , rjM) = -(sin^,r/^(^)). 

Similarly, Cr] = if is odd. 

We denote the even and odd subspaces of E by Ec and Eg respectively, to remind us 
that the even (odd) tangent vectors correspond to cosine (sine) perturbations of the order 
parameter. The even-odd decomposition also respects the locked-drifting decomposition, 

so we can express E^ = E\.® E\ and E'^ = E^ © E^ as the direct sum of even and odd 
subspaces. Since any even or odd 77 e is completely determined by rj^ for a; > 0, we can 
identify 

E^^ El^ L\n^_^_,C) and E^ ^ E^ ^ L^^l^W) 

where = fl [0, 00), Q'i^ = Q'^ n [0, cxo). The operators M and L both preserve the 
even and odd subspaces. (Proof: if i] is even, then v^rji^j is odd, so D{v^r]i^) is even. And 
[sin9)pi_j{9) is odd, so D {{sin 9) pi_j{9)) is even; hence Mi] and Lr] are even. The proof for 
T] odd is similar.) Hence we can analyze the spectrum of L restricted to the subspaces Ec 
and Eg separately, and combine the results to obtain the spectrum of L on 
M has only continuous spectrum, and it's fairly easy to describe: 
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Proposition 3. a{M) = [-KR, -y/{KR)^ - 1] if p is fully locked, and a{M) = [-KR, 0]U 
Rz if p is partially locked (here IRz denotes the imaginary axis). All A e cr(M) are in tie 
continuous spectrum. 

The proof of this proposition is somewhat technical, so we postpone it to the end of 
this section so as not to interrupt the main thread we are developing. When two closed 
operators on a Banach space differ by a bounded operator of finite rank, as is the case for 
L and M, their spectra are closely related. In fact. Proposition 3 is the key ingredient in 
establishing the main result of this paper. We present the proof of this theorem below, 
with some of the steps deferred to Sections 7 and 8. 

Theorem. The spectrum (j(L) consists of cr{M), an eigenvalue at 0, and perhaps one other 
eigenvalue A e [— •>/ {KRY — 1, 0) if p is fully locked. Except for and \, the spectrum of 
L is all continuous. 

Proof. Suppose A G a"(M). The spectrum of M is all continuous, so the image of the 
operator XI — M, which coincides with the domain of the unbounded operator {XI — M)~^ , 
must be dense in E and have infinite codimension (otherwise we could extend {XI — M)~^ 
to a closed operator defined on all of E by adding a bounded, finite rank operator, but 
then the closed graph theorem would imply that (AI — M)~^ is bounded on E). Therefore 
AI — L is also not onto, so A e cr(L); hence a{M) C cr(L). Reversing the roles of L and 
M in this argument shows that ac{L) C a{M). So to complete the proof we need to 
prove three things: is always an eigenvalue for L (this is not surprising considering the 
rotational symmetry of the Kuramoto model); L has at most one other eigenvalue A 7^ 0, 
which satisfies A G [—a/ {KRY — 1, 0) and occurs only in the fully locked case; and the 
residual spectrum (7r{L) = in all cases. 

Let's begin with the eigenvalues. As discussed above, it suffices to study L separately 
on the even and odd subspaces Ec and Eg', on Ec, L is given by 

{Lt])^ = M^T]^ + K{Cr])D{{sine)p^). 

M is defined pointwise as a function of a; G and as we shall see in the course of the 
proof of Proposition 3, for any given A G C the operator {XI — M^) is invertible for all 
but at most countably many values of uj. So suppose r] E Ec is an eigenvector for A. We 
must have Cr] 7^ 0, since otherwise A would be an eigenvalue for M. So we can assume 
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that Crj = K ^ if we like. Then rj is determined uniquely by the formula 

Conversely if A ^ a{M) then this equation defines an element i] G Ec- But if A G a{M) 
then r] may not be in Ec, since i] may fail to be integrablc w.r.t. the density function g{u). 
X is an eigenvalue for L on Ec if and only if rj E Ec and A satisfies the self-consistency 
relation Crj = K~^. Define the function he by the formula 

hc{X) = [ {cos e, {XI -M^y^D {{sin e)p^))g{io)dio; 
Jo. 

the domain of he is defined to be those A G C for which the integrand above is integrable. 
Then a necessary condition for A to be an eigenvalue for L on E^ is that A satisfies the self- 
consistency equation hc{X) = K~^; this condition is necessary and sufficient if A ^ a{M). 
The situation is similar for Es. Here L is given by 

(L?7)^ = M^?7^ - K{S'n)D{{cose)p^). 

We define 

hs{X) = - [ {sine,{XI - M^)-^D {{cose) p^))g{u)duj; 
■In 

if A is an eigenvalue for L on Eg, then A must satisfy the self-consistency equation /is (A) = 
K~^. We call he and hg the characteristic functions for L. We will derive explicit formulas 
for these functions in the next section, and then in Section 8 we shall prove 

Proposition 4. The equation hc{X) = has at most one nonzero root X, and only 

if K > Ki. This root satisfies — < A < 0. In addition, X = is a root if 

and only if K = K^- The equation hs{X) = has A = as its only root in all cases. 
Furthermore, the roots of the characteristic equations are in fact eigenvalues of L on E^ 
and Eg respectively. 

This completes the description of CFp{L). Notice that we have to be a little careful here: 
a root A G a{M) of one of the characteristic equations is not automatically an eigenvalue 
of L, since the associated eigenvector ry might not be integrable w.r.t. g{oj). Fortunately 
as we shall see in Section 8, this doesn't happen. 

Now what about the residual spectrum of L? To answer this, we'll need a concrete 
description of the elements v G L^(f2, VF)*. Let Wq = L'^{S^,C) n 1^ be the subspace 
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of W consisting of measures of the form (f)dm, where is on S^. The subspace Wq 
naturally has the structure of a separable Hilbert space, and its image in W is dense, which 
implies that the natural map L^{n,W)* L^{n,Wo)* is injective. Now L^(n, VFq)* = 
L°°(0, Wq) (see Lang [1993, p. 188] for this result), and the Hilbert space Wq can be 
identified with the space of functions on with integral 0. So an element u E E* can 
be represented as a function c<; i-^ z/^^, where is a function on for each frequency uj. 
If ry e the pairing ((77, v)) of ry and v is given by integrating the function uj 1-^ (I'ojjfjui)- 



Jn 



Now suppose the operator XI — L does not have dense image on Ec- Then there must be 
a nonzero element v E E* such that 



(((A/ - L)r7, u)) = {{{XI - M)r^, v)) - K{Crj){{D{{sm e)p) ,u)) = 

for all 1] E Ec in the domain of M. Since M has only continuous spectrum, we cannot have 
((D((sin^)p), jy)) = 0, so we may assume that ((D((sin^)p), i^)) = if we like. Then 
we have 

(((A/-M)r/,i.)) = Cr/ 
for all r] E Ec in the domain of M, or equivalently, 

{{ri,u)) = C{XI-M)-'rj 

for all r] E Ec in the range of XI — M, which is dense in Ec- The operator (A/ — M^)~^ 
exists and is bounded for almost all u>, and so we have 

/ {iyM,vM)9i^)duj= [ {cose,iXI - M^)-'riUO))gic^)du 
Jo, Jq 

= [ {{XI-M*)-^cose,rjU0))9i^)dw 

for all r] E Ec in the range of XI — M, where M* denotes the adjoint of the operator M^j. 
This uniquely determines i/: 



u^{d) = {XI -M*)-'^ cose 
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for almost all uj. Now 

{{D{{suid)p),v))= [ {{XI-M::,)-^cosd,D{{sme)p^))g{uj)duj 
Jo. 

= I {cos9,{XI -M^y^D {{sin 9) p^))g{co)dw = hc{X) , 
Jq 

so the self-consistency equatfon ((D((sin ^)p), i/)) = is exactly the same as before: 

hc{X) = K~^. The same argument applies to L on Es. But by Proposition 4, any roots of 
the characteristic equations are eigenvalues of L, so the only way XI — L can fail to have 
dense image is if A e crp{L), and thus we conclude that ar{L) = 0. 

QED 

We now turn to the proof of Proposition 3. 

Proof of Proposition 3. We'll analyze the spectrum of M on the spaces and E'^ 
separately and then combine the results. The most general tangent vector r] ^ E^ has the 
form 

Vuj = c{uj)DS0^ 

where the coefficient c is an L^-function of a; G w.r.t. the measure g{u)du. Let (f) be 
any smooth function on S"^. Then 

(0,M,r7,) = {<p, -D{v^r]^)) = {v^cj)' ,r]^) = -c{u)M)\e^). 

But Vt^{9tS} = 0, so 

{(f>,M^r]^) = -c{uj)vU0u>)(l>'{9^) = c{uj)KRcos9^cf>'{9^), 

and hence 

Muj-rjoj = -KRc{oj) cos 9^D5e^ ■ 

In other words, M is just multiplication by the function —KRcosO^. This explicit de- 
scription shows that M is a bounded operator on i?', since the function cos^tj is bounded. 
{XI — M^) is multiplication by (A -|- KRcos9^) which is nonzero a.e.; hence XI — M has 
trivial kernel on EK 

The inverse oi XI — M on E^ is multiplication by (A -|- KR cos 9 this operator 
is bounded on E^ = L^(1]',C) if and only if the function a; i-^ (A -|- KRcos9^^)~^ is 
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essentially bounded. This is equivalent to A not in the image of the continuous function 
uj I— > —KRcosdoj as uj ranges over the locked oscillators. Hence the spectrum of M on 
is [-KR, -^{KRY - 1] if p is fully locked, and [-KR, 0] if p is partially locked. For 
these A the operator (A/ — M)^^ is unbounded, but is defined on coefficient functions c 
that are supported away from the roots of A + KR cos 9^, and hence is densely defined. 
Therefore the spectrum of M on E'' is all continuous. 

When the state p is partially locked, we will show that the spectrum of M on E^ is Wi 
and is all continuous, and this will complete the proof. For oj e Cl^, the operator is the 
adjoint of the first-order differential operator on defined by N^cf) = v^cf)'. Note that 
is a closed, unbounded operator on the space of complex- valued functions C''{S^,C), 
and is a bounded operator from C'^{S^,C) to C''~^{S^,C). We will see that for each oj 
with |a;| > KR, the spectrum a{N^) = ap{N^) = {2TvinCuj \ n G Z}, where Ceo is the 
coefficient given by (8). So if Re A 7^ 0, then {XI — N^)~^ exists and is a bounded (in fact 
compact) operator on C''{S^, C). Therefore (A/ - M^)'^ = ((A/ - N^)~^)* is a bounded 



operator on C^{S^X)*, with image in C'=-i(5\ C)*. Notice also that {XI - N^)l = A, so 
{XI-N^)-^1 = X-^. Hence ii rj^ e C''{S^,Cy satisfies {l,r]^) = 0, then 



Therefore {XI - M^)'^ preserves the subspace W C C^{S^,C)\ since C'=-i(5\C)* C 



The operator N^^ depends continuously on a;, and hence so does the operator {XI — 
Mj^)~^. We will also show that the norm of {XI—N^j)~^ (as an operator mapping C^{S^, C) 



this proves that if Re A 7^ 0, then XI — M has a bounded inverse on E'^ defined for r] e E' 
by the rule 



If A e K but A 7^ 0, then {XI — M^) ^ is defined at all frequencies except ±0;^, where 



V 

is obtained by solving A = 27riCu,n for Un- Note that the problem frequencies ±0;^ have 
no limit point in the open intervals ±{KR, 00). {XI — M)~^r] is defined for any r] E E^ 



(1, {XI - M^)-i ry^) = {{XI - N^T^r^^) = X-\l, ry^) = 0. 



to itself) is bounded as a function of u. Since 



{{xi-N^)-^y 



||(A/-iV,)-i||, all 



{{XI - M)-^r,)^= {{XI - N^)-y 
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that vanishes on a neighborhood of {±a;„}, and the set of such rj is dense in E^. This 
proves that {XI — M)~^ is defined on a dense subspace of E^; in other words, the image 
of A7 — M is dense in E^. However, we claim that A7 — M is not surjective, and hence 
A is in the continuous spectrum of M on E'^. To see this, let ^ be a non-trivial solution 
to {XI — N^^)(f) = for some problem frequency uj^ {cp is guaranteed to be C°°). Choose 
rj e E'^ such that ((/>, rj^) = 1 for all uj in some neighborhood of LVn (we can do this because 
(f) is not constant), and suppose {XI — M)~^r] exists in E^. Then 

{XI - N^)^ = {{XI -N^J + {N^„ - N^)) ct>={un- u;)<P', 

and so for uj sufficiently close to a;„ we have 

1 = (0, Vo.) = ((A/ - N^)cf>, {XI - M^y'fj^) = {LOn - u;){ct>', {XI - M^y'rj^). 

Therefore 

1< |a;-a'n|||</''||||(A/-M^)-Vll 

as w — > a;„, where the norms of 0' and {XI — M^)~^r]^ are taken in the spaces C'^{S^,C), 
C''{S^,C)* respectively. But this implies that the function oj i— > ||(A/ — M^)~^r]^\\ is not 
integrable, which is a contradiction. 

Special care must be taken when A = 0, since e crp{N^^) for all cj. The spectrum 
of a closed operator is always a closed set, so is in the spectrum of M on M has 
trivial kernel in E^, since the only solutions to the equation D{vuj'r]^) = are multiples 
of pc^, but any i] E E'^ satisfies {l,r]i_j) — for all to. We can invert M,^ on the subspace 
W C C"=(5^C)* as follows: 

Vu; = -v^ [D r]^ —J ^ — m , 

V {vui ,m} J 

where D~^rj^ denotes the unique distributional antiderivative of rj^ on determined by 
the requirement {l,D~^r]^) = 0. If 77^; = for uj in some neighborhood of rkKR, then 
M~^r] is perfectly well-behaved, so is defined on a dense subspace of E*^; in other 

words, the image of M is dense. Therefore we see that M has empty point and residual 
spectra on i?^, and continuous spectrum Mi, as desired. 

Now let's complete the analysis of the operator N^. Observe that A is an eigenvalue 
of if and only if the equation 

X(f) — Vu,(f)' = 
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has a non-trivial solution on S^. This equation has general solution 



0(^) = A; exp (^xj^ v^ry^dr^ , 



and assuming the constant k ^ this function is periodic if and only if exp(^) = 1, which 
is equivalent to A = 27riC^n for some n G Z. Therefore (7p[N^^) — 21^10 JIj. If A ^ ^TiiCJIj 
then we rewrite the ODE 

\^ - v^iP' = (t> (15) 

as 

We can solve this equation with the aid of the integrating factor 



= exp ^-A J v^{t) ^dr^ ; 



the unique periodic solution is 

m = -iM-' (^j\M)vUr)-'Hr)dT + {^U'^7r) - l)'' J%Ur)vUr)-' ^(r) dr^ . 

7t^(27r) = exp(— ^) 1, so this function is well-defined. The right-hand side defines the 
bounded operator {\I - N^)'^ from C''{S\ C) to C''+\S\C), and so a{N^) ^ (Jp{N^) = 
{27TinC^\ne Z}. 

Next we establish the necessary norm bounds for {XI — iV^)~^. The differential equa- 
tion (15) has an irregular singular point at the limiting values u) = ±KR, since the function 
= uj — KRsinO has a double root at ^ when u = KR (and a double root at— ^ when 
uj = —KR). So the fact that ||(A/ — A^a;)~-^|| is bounded as a function of uj is not trivial. It 
suffices to prove the existence of constants C/., A; > 0, depending only on A and KR., such 
that 

||(A/-iV^)-VlU <Cfc||(/>||fc 

for any e C''(5'^C), and any oj with |a;| > KR. For A; = 0, we take e C(S'\C) and 
put ij) = {XI — N^)~^(j) e C-^(S'-'^, C). Suppose |V'(^)P has a positive maximum at ^ = ^o- 

^Q\m? = i^\e)m+mno) = 2Re{i;'{e)W)i 
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and so Re(V''(^o)V'(^o) = 0. Multiply (15) by V'(^o) and take real parts of both sides to 
obtain 

(ReA)|V'(^o)|' = Re((/>(^o)VW); 

hence 

\^|;{^o)\ < |ReA|-V(^o)|, 

and so we have HV'llo < (^0110110 with Cq = iReAl"-*^. This proves the base case /c = of our 
assertion. 

Now we proceed by induction, and assume the existence of the constants Cq, Ci, . . . , 
has been established. As before, let ip = {XI — N^)~^(f)^ where now e C^'^^{S^ and 
hence V' £ C'^+^(S'^, C). Differentiating (15) k times gives 

The terms v'^~^^ are all bounded by KR^ so by induction the right-hand side is bounded 
by C"||(/)||/;, and hence C"||0||fc_|_i, for some constant C . If I^/jI'^+^II takes its maximum at 
9 — 9o, then we obtain 

U^''+%<C'\uj-KRsmeo\-^mk+i. 

Now we differentiate (15) one more time and imitate the proof for the case k = 0. We 
write the result in the form 

(A - (A: + l)v'Ji;y''+'^ - = 0^'+'^ ^ ^ d ^)^- 

and observe that the right-hand side is bounded by C"||0||fe+i for some constant C". 
Multiply by ■^['^+^l(^o) and take real parts of both sides to obtain 

||V^['^+il||o <C"|ReA+(A;+l)Mcos^o^'||0lU+l• 
Fortunately, the function 

m{d,uj) = min (C"|a; - Msin^|-\ C"|ReA + (/c + l)Mcos^|-^) 
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is bounded for 9 e S^^ \uj\ > KR. (If not, we could construct a sequence (a;^,^^) such 
that m{uj„, $n) — > oo, but this imphes that both 

I sin6'n| 1 and cos6'„ -7rT~TTW^ ^ 0, 

[k + 1)KH 

which is impossible.) Therefore we have | IV'^'^"'""^^ ||o < (^"'II^IU+i, where C" is any bound 
on m{9, Lo), so we can take Ck+i = C" + Ck, and we're done. 

QED 

Remark on cr(L) for non-special positive states. In this case, we can carry out a 
similar analysis of the continuous spectrum (Tc(M), but now it will contain positive real 
values, coming from the essential range of the function uj i— > —KRcosO^ = KR cos 6,^. 
The same reasoning used in our main theorem shows that these values are also contained 
in (t(L), so these states are not linearly stable. We omit the details since these states are 
of minor importance to us. 



7. Characteristic Functions for L 



Our next task is to compute the characteristic functions he and hg-, given by 
hc{X)= I {cos9,{XI -M^y^D {{sine) p^))g{oj)doj, 

hs{X) = - I {sine, {XI - M^)-^D {{cose) p^))g{io)dw. 
Jo. 

To do this, we'll have to compute {XI - M^)-^D {{sine)) and {XI - M^)-'^ D {{cos e) ex- 
plicitly. Let's begin with the locked frequencies cu e fl''. As we saw in the previous section, 
for these frequencies {XI — M^^)"^ is just multiplication by (A + Ki? cos 6*^; ) ~ ^ . Hence we 
obtain the contributions h^. and h[ resp. from the locked oscillators to the characteristic 
functions he and hs'- 

h[{X) = [ {X + KRcose^)-^{cose,D{{sine)p^)) g{u;)dw 
= / {X + KRcoseu,)~^{sine,{sine)pu,)) g{uj)dw 



sin^ e,. 



■ g{uj)dw 



(16) 



/q; a + KR COS e^j 

r-KR , ,2 
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and 

/i^(A) = - / {X + KRcos9^)-^{sme,D{{cose)p^)) g{u)div 

g{uj)(hj (17) 



cos^ 6^ 



Qi A + KR cos 9^ 

{KRf - 



(KR)' J- 



: g{uj)dui). 



-KR A + y(M)237:;2 

The functions /ic(A), h[{X) are (up to the constant K) the N ^ oo hmits of the rational 
functions i?s(A), -Rc(A) resp. that we defined in our study of the finite- Kuramoto model 
[MiroUo and Strogatz 2005], as was to be expected. After all, completely locked states are 
just the infinite- A'^ analogues of fixed points for the finite- A'^ Kuramoto model. However the 
partially locked states have no finite-AT analogues, so we should expect to see something 
new there. 

Now suppose p is partially locked, and consider the drifting frequencies oj E 0,^. For 
almost all uj e Q^, the distribution {XI — M^)~^D {{sin 9) p^^ will be a smooth measure of 
the form ai_j{9)d9, where a^j is a smooth function on S^. The operator A7 — M^^ applied 
to the measure a^(9)d9 gives the measure (Xaa,(9) + (vrjarj)' (9)) d9 , and D((sm9)pi^) is 
ju.t the — {^y^. The.efo.e must .atMy the fl.st-oMe. ODE 



XaUO) + {vMaM)' = 



sin^ 

MO) 



Our strategy for solving this equation is to express a^, in the form P^/v^, so as to clear 
out the denominator above. The corresponding equation for is 

XpU9) + v^{9)PM-vMPM = CUv^{9)cos9-vU0)sm9). (18) 

Now a fortunate miracle occurs. The functions sin^, cos^ and v^{9) are all contained 
in the vector space V spanned by 1, cos^ and sin^. Notice that V is closed under the 
operation 

so we restrict our search for solutions to (18) to functions of the form 

Pui{0) = co{u!) + ci(a;) cos ^ + 02(0;) sin ^, 
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with coefficients Cj that depend on uj. (The miracle is that the three-dimensional subspace 
consisting of measures pijj of the form (^Pi^j /v^^dO , with Pijj being a linear combination of 
1, sin^ and cos^, is invariant under the differential operator M^.) Substituting Vui{0) = 
u) — KRsinO in (18) gives 

{Xco{uj) + KRci{uj)) + {KRcq{uj) + Aci(a;) + ujc2{oj)) cosd + {—lvci{lv) + Ac2(a;)) sin^ 
= (jjC^ cos 6. 

We equate coefficients and solve for Ci{u>) to obtain 

KRwa 



co(cj) = 
ci{u!) = 

C2(<^) 



X^+uj^-iKR)^' 
A2 + a;2 - (Ki?)2' 



A2 + - {KR)^ 

So the contribution to he from the drifting frequencies is given by 
^c(A)= / (cos^,(A/-M^)-ii^((sin^)p^))^(a;)cL; 

cos6'(co(w) + ci{u>) cos 6* + C2(w) sin6')f ^^(6*) d9 ) g{u>)dw 

C0S2^ 



\uj\>KR 







f . , f f'^ cos' 9 
= / Ciiuj) / 7 — — T^TTT do q{u)du} 



\w\>KR 

(the other two integrals vanish because the integrands have periodic antiderivatives on S'-^). 
We can evaluate the inner integral using integration by parts: 

cos^ „ 1 /•^'^ sine 

df) 



cose n.nj_r 

Jo (a;-Msine)2 KRJo u - KRsind 



{KRY Jo \uj-KRsine 



i\de 
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Therefore 



\u)\>KR 

= WW I A^+a;^-(i^i?)^ ^^'~'^^^-^^^^^^ (19) 

\uj\>KR 

\uj\>KR 

We repeat this procedure to compute the smooth measures (A/ — Mi^)~-'^D((cos 9)pS) 
for uj e O'^; express this measure in the form j3^{9)d9, where now (3^ satisfies the equation 



We put f3u>{Q) = Co (a;) + ci(a;) cos^ + 02(00) sinO and solve for the coefficients Cj as before, 
this time obtaining 

, , KRXC^ 
Co (a;) 



ci (uj) = 

C2{uj) = 



C^{oo^-{KRf) 
A2 + a;2 - [KRf ' 

Aa;Ct^ 



A2+a;2-(M)2- 
The contribution /i^ to hg from the drifting frequencies is given by 



hi{X) = - I {smd,iXI - M^)-^D{{cosd)p^))g{uj)dw 

= - J sin6'(co(a;) + ci(a;) cos^ + C2(a;) sin6')t;J^(6') (i6'^ 5f(a;)(ia; 

\uj\>KR 

J Co (w) sin 6* + C2(cc') sin 9 



{u-KRsvney 

\lj\>KR 
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Substituting the coefficients cq and C2 gives 

^ar^, f AC^ f f'^'' -KRsm6 + usm''6 \ ..^ 

^'(^)= J X^+.^-iKRr [Jo {co-KRsinOr 'T^"^'" 

\u}\>KR 

- J X^ + u-^-iKRY U U - Msin^ " (a; - Msin^)^ ) ^^""^ 

|a;|>i<rit 



(20) 

So the characteristic functions he and /i^ are given by hc{X) = /^[.(A) + h'^{X) and /is (A) = 
^i(A) + /if (A), with the functions h^, hf, h\ and given by (16), (17), (19) and (20) 
above. This completes the derivation of the characteristic functions. 



8. Roots of the Characteristic Equations 

Our final task is to prove Proposition 4, which we restate here for convenience: 

Proposition 4. The equation hc{X) = has at most one nonzero root X, and only 
if K > Ki. This root satishes — ^/(K^p"^^ < A < 0. In addition, X = is a root if 
and only if K = K^- The equation hs{X) = has A = as its only root in all cases. 

Furthermore, the roots of the characteristic equations are in fact eigenvalues of L on Ec 
and Eg respectively. 

Proof. Notice first that 

and /if (0) = in the partially locked case, so A = is always a root of the characteristic 
equation /is (A) = K~^. This is no surprise, since the rotational symmetry of the Kuramoto 
model implies that A = is always an eigenvalue of L. The corresponding eigenvector is 
just Dp e Eg-, we have S{Dp) = —R, and so 

L{Dp)^ = -D{v^Dp^ - KR{cose)p^) 

= -D\v^p^) = 0. 
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Next, let's consider the fully locked case, so ft = [—1, 1], K > Ki (equivalently, KR > 
1) , he = h^, and hs = hg. If A e C has nonzero imaginary part, then the same is true for 
hc{\) and hs{X); to see this, multiply the numerator and denominator in the integrands for 



/i^(A) and hg{X) by A+a/ (KR)'^ — u"^. Hence the characteristic equations can have only real 
roots in this case. So let A G M. If A < —KR then hdX), hs{X) < 0, so A is not a root of the 
characteristic equations. (Note that hc{—KR) is well-defined, whereas hs{—KR) — —oo.) 
hc{X) and /is (A) are undefined for -KR < A < -^{KRf - 1, since the integrands in the 



formulas for hc{X) and /is (A) have simple poles at a; = iy^ (KR)'^ — A^ G (—1, 1). The 



functions he and hg are both defined and positive on (—^/{KR)^'^^, oo), but might take 
the value +00 at A = —■^/{KRY — 1- And both functions are strictly decreasing on the 
interval [— — l,oo), so the characteristic equations each can have at most one 
root here. Since we already saw that /is(0) = K~^ in all cases, we conclude that there are 
no other roots to the characteristic equation hs{X) = K~^ in the fully locked case. 

Next wc claim that hc{0) < K~^, which implies that any root A of the equation 
hc{X) — K~^ must satisfy A < 0. To see this, observe that 

KR ( J, 



_ 72 



2 ' 



The function h{s) = — \/l — changes sign from negative to positive at sq = 

and g is non-increasing on [0, 1] by assumption. Therefore 

J (/ic(0) - i^-^) = / h{s)g{KRs)ds+ [ h{s)g{KRs)ds 

2 Jo Jso 

<g{KRso) / h{s)ds + g{KRso) / h{s)ds 

Jo Jso 

= g{KRso) ( yj=f - Vl - ds = 0. 

This inequality is strict except in one special case: when KR = 1 and g is constant on 
[0, 1] ; recall that in this case the two critical coupling constants K^. and Ki are equal. So 
in the fully locked case with K > Ki the equation hc{X) = K~^ has at most one root, 
which satisfies — {KR^ — 1 < A < 0. li K = Ki then the equation hc{X) = K~^ has 
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no roots, except for A = in the special case where g is constant on [0, 1], or equivalently 
when K = Ki = Kc . 

li K > Ki, then the characteristic equation hc{X) = has a root in the interval 



[-^{KR)'^ - 1, 0) if and only if K{-sJ{KRY - 1) > R-^ . If we examine the formula for 
hc{ — \J (KR)'^ — 1), we see that we can make this value as small or as large as we like by 
varying the density function (to make hc{ — \/ (KR)"^ — 1) small, take the density function 
to be concentrated at w = 0, where the integrand vanishes). So it is possible that he has 
no roots in [-y/{KR)^ - 1, 0). If the endpoint A = -^{KR)'^ - 1 e a{M) happens to be 
a root, then we need to check that the associated eigenvector r] given by 

V. = ^{^+V{KRr-^'r'D5o^ 

is in fact an element of Ec- The coefficient function above becomes infinite as c<; — > 
±1, but is nevertheless integrable against g{u>), since the integrand in the formula for 
hc{—-\/ {KRY — 1) has the same asymptotic behavior as a; — > ±1. So is a bona fide 
element of Ec-, and A is an eigenvalue of L on Ec- The same thing happens in the special 
case with KR = 1 and g constant on [—1, 1]: we have 

Vl — u;^ 

and the coefficient function is integrable on [—1, 1], so A = is an eigenvalue of L on E^.. 
This completes the proof for the fully locked case. 

Now we turn to the partially locked case K < Ki, and prove that both characteristic 
equations have no roots A 0. Let's begin with A e R. If A < —KR, then hc{X), 
/is (A) < 0, and if —KR < A < 0, then hc{X) and /is (A) are undefined. To complete the 
proof for A e R, we wiU establish that hc{X) < K~^ and hs{X) < K~^ for all A > 0. To 
simplify the notation a bit, replace A with KRX; then the first inequality is equivalent to 

/ (- ^== - Vl - sA g{KRs)ds + r ,„ \ , ^=== g{KRs)ds < 

for all A > 0. As before, the function h{s) = ^^^^_^2 ~ \/l — -s^ changes sign from negative 
to positive at a unique point sq e (0, 1), so 
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the integral on the right is negative for A > and g is non-increasing, so we have 

For the second integral, observe that 

/ T^^ o ; g(KRs)ds < g(KR) / — ; ds 

(strict inequality holds because g{KR) > and g{(ju) — > as a; — > oo). So it suffices to 
prove that 

( y/l^] ds + /"t^^-A ^ cig<0 

Jo VA + yr^ J Jl A2 + s2_l ^ + ^^23T - 

for all A > 0. These integrals can be evaluated explicitly, and in fact equality holds above: 



/ , ds + [ —: — ^ ■ ^ ds — f \/ 

Jo X + VT^ Jl A2 + S2-1 s + V^^ Jo 



1 — s'^ ds = — 
4 



if Re A > 0. The relevant explicit formulas are 



and 



/ g ds^--X^ + X+^ + A VA^ - 1 tan-^ y/X^ - 1 

Jo A + Vl - 2 4 

/•oo \ 

/ -^-4^ y ds = lx''-X- Ay^A^tan-^ a/a^ 



for A > 1. Keep in mind that the integrals above are analytic functions on the domain 
Re A > 0, so proving they agree for A > 1 is sufficient. 

The proof of the inequality hs{KRX) < is similar: we need 



/ == - v^r^ g{KRs)ds + / , , • "^'r^ g{KRs)ds < 

Jo VA + VI - s2 J Jl X^ + s^-1 s + Vs^ -1 



for all A > 0. We have 



/ ( ^ , - ^/l-sA g(KRs)ds = -A / g(KRs)ds 

Jo \X + VT^ J Jo X + VT^ 

< -Xg(KR) [ ~ ds 

- 'Jo A+vr^ 
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and 



as before. This time we find that 



poo V 1 



TT. . A 



2 

ds — I — - — , ds 



2 



= --A + 1 + , tan-i VA2 - 1 

2 VA^^ 

for all A > 1, and this completes the proof that the characteristic equations have no 
non-zero real roots. 



Now let's rule out any complex roots. Since hc{X) = hc{X) and hs{X) = /is(A), it 
suffices to prove that ImA > implies Im/ic(A) < and Imhs{X) < 0. Let ImA > 0. If A 
is pure imaginary, then hc{X) and /is (A) are defined only if O = [—1, 1] and ImA > KR, 
and it's easy to see that Im/ic(A) < and lmhs{X) < in this case (all the integrands 
have negative imaginary part). So assume Re A 7^ 0. Replacing A with KRX as before, we 
see that it suffices to prove 

/■I §2 r°° X s 
Im / ^ — g{KRs)ds + Im / ^ ■ — — ^ ^ g{KRs)ds < 

and 



A + ^r^ ' ' 7i A2 + s2_l ^ + ^^2^ 



for ImA > 0. The proofs of these two inequalities are identical, so we'll present only the 
first one. ^ 

Im ( ^ ) < 

VA + ^T^y 

for < s < 1, so 

/ 'I 2 9{KRs)ds < g{KR) Im / ^! ^ ds. 

The function 

Im ( 7i ^ = ImA ^ ^ ^ 



A2 + s2 _ 1 y \^ 1^2 _^ 
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changes sign from negative to positive at sq = a/|A|M-T. So by the same argument as 
before, 

A s f°° A s 
Im / ;s , q(KRs)ds < q(KRso) Im / , ds. 

Therefore 

/ X , r, 2 9{KRs)ds + Im / ■ ^— ^ g{KRs)ds 

Jo A + Vl - Ji A^ + - 1 s + - 1 

< g{KRso) Im / , _ rfs + g{KRso) Im / ^ ^ r , „ ^ 

/ /■! g2 f°° X S \ 

The right hand side is if Re A > 0, so we're done in this case. If Re A < 0, let A* = —A 
be the reflection of A in the imaginary axis. Then 

Im I , 1 < Im 



X + VT^J \X* + VT^ 

for < s < 1 and 

Im ( — I = Im 



X^ + s^-lJ V(A*)^ + s2-l, 

for all s > 1, so we are done by the previous case. 

To finish up, we need to see when A = is a root of the characteristic equation 
hc{X) = K~^. The same argument as in the fully locked case shows that this holds if and 
only if g is constant on [—KR,KR\, which is equivalent to = Kc- Since G a"(M), 
we need to check that the associated eigenvector r] given by the formulas in Section 7 is a 
bona fide element of the tangent space E. If |a;| < KR, rj is given by 

and if |a;| > KR, rji^{0) = a^{9)d0 with 

1 |a;| — Ki? + a;sin^ 



27r ya;2 - (KR)^ - KRsmd)^' 
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The function to i— > ((i^i?)^ — a;^)~^ has an integrable singularity at a; = ±KR, and 
fortunately, the same thing happens for the drifting frequencies. To see this, express 
uj = KR{1 + e), with < e < 1, to analyze the singularity at a; = KR (the singularity at 
—KR has the same behavior). Then 

1 1 + f -l + (l + e)sin6' 1 1 + ^ ^ . ^^^/' 

"-^'^ = 2^ • TlTTl^ ■ (l + .-sin^)^ =2^ -717^1^^^°^^' + '"""'^^ • 

If (/) is any C°° function on , then 



1 1 + e '•^'^ 



Hence Hry^jH, taken in C'^{S^)*, satisfies 



/ 0"(6')log(l + e-sin 6')d6' 
Jo 



c Z-^" 

lha.ll < ^ / |log(l + e-sin6')|d6' 
Jo 

for some constant c, and the integral above is bounded as e — (this follows from the fact 
that log a; has an integrable singularity at 0). Hence 77 is an element of E, and the proof is 
complete. 

QED 

The eigenvector described above has an interesting interpretation. As we saw in 
Section 4, if the density function g is constant on [0,a;o], the model with critical coupling 
Kc = 2/'Kg{Qi) has a family p{t) of special positive states parametrized by t e (0,a;o]. The 
tangent vector 77 described above is exactly the derivative of p{t) with respect to t, at the 
value t = KR. 

Remark on cr(L) for the incoherent state. All our methods apply as well to the 
incoherent state, and are in fact much easier to work through in this case. The details are 
in [Strogatz and MiroUo 1991], so we present here a brief summary of the results in the 
context of our current formulation. There are no locked frequencies now, and the operator 
is just —ujD on S^, which has spectrum iuj'L. So M has purely continuous spectrum 
Mi, and the characteristic functions he and hs are given by 

1 f°° A 
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The only possible roots of the characteristic equation are real and positive, so the contin- 
uous spectrum of L is also Rz. If A > 0, we express 



SO the characteristic equation hc{X) = has a root A > if and only if > Kc, and 
hence the incoherent state is linearly unstable in this case. It is of course no coincidence 
that the incoherent state loses stability at the same critical coupling at which the positive 
fixed states are born. 

9. Concluding Remarks 

We can now completely describe the spectrum of L for all special positive states in the 
infinite- AT Kuramoto model (Figure 2). For the fully locked case with K > Ki the spectrum 
consists of the closed interval [—KR, —^/{KR)'^ — 1], at most one negative eigenvalue in 
[— ^ (KR)"^ — 1,0), and the eigenvalue at coming from the rotational symmetry of the 
Kuramoto model. This result is not surprising; the corresponding finite- A*" model [MiroUo 
and Strogatz 2005] has stable fixed points for most choices of frequencies uji as N —>■ oo, 
provided the coupling K > Ki. The eigenvalues for these fixed points are all negative, and 
so cr(L) is in some sense just the limit of the spectrum in the finite- A?" case. 

When we make the transition to the partially locked case at K < Ki, the spectrum 
explodes to include the entire imaginary axis, along with the segment [—KR, 0], as shown 
in Figure 2c. The presence of negative values in cr(L) suggests some sort of asymptotic 
stability, but the presence of Rz in the spectrum suggests more neutral behavior. So in 
some sense what we have accomplished is to rule out any kind of exponential convergence 
to fixed states in the partially locked case. 

On the other hand, the order parameter R{t) might still be able to approach its 
stationary value exponentially fast, due to a phase- mixing mechanism akin to Landau 
damping; this is known to occur for the incoherent state under some conditions [Strogatz et 




to see that he is strictly decreasing on (0,oo). Furthermore 
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al. 1992; Balmforth and Sassi 2000]. One should also keep in mind that the partially locked 
states have only been shown to be linearly neutrally stable. The small nonlinear terms 
neglected here could therefore prove crucial to understanding the full stability properties 
of the fixed states. A careful nonlinear analysis along these lines, ideally one that is global 
in character, may be the next logical step in the ongoing attempt to make sense out of 
Kuramoto's marvelous calculation of 1975. 
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Figures 



Figure 1: Order parameter R for fixed states of the infinite- A?" Kuramoto model, as a 
function of the couphng strength K. Points {K, R) on the curve C correspond to special 
positive fixed states. Points {K, R) below the curve and above the K-axis correspond to 
positive fixed states with nonzero weight functions w; all of these non-special states turn 
out to be unstable, as we remark at the end of Section 6. The curve C has a vertical 
segment if and only if the density function g is locally constant at 0. And when O, the 
support of g, is [—1,1], points on or above the hyperbola KR — 1 correspond to fully 
locked states, and points below correspond to partially locked states. We present three 
qualitatively different cases: (a) Q = [—1,1] and g{x) < g{0) for x ^ 0. (h) fl = [—1,1] 
and g is locally constant about its maximum. Here C has a vertical segment at K = Kc, 
corresponding to a 1-parameter family of special positive states, all with K = Kc but 
different values of R. (c) = R and g{x) < g{0) for x ^ 0. Full locking is never achieved. 

Figure 2: The spectrum a{L) for the special positive fixed states, (a) For fully locked 
states with K > Ki, X and are eigenvalues. The rest of cr(L) is continuous. The zero 
eigenvalue follows from the rotational symmetry of the Kuramoto model. In contrast, the 
eigenvalue A is not present in all cases; it exists if and only if hc{-^/{KRY - 1) > K'^, 
as shown in the proof of Proposition 4. (b) For fully locked states at the bifurcation value 
K = Ki, the spectrum contains an eigenvalue at 0, and the rest of cr{L) is continuous, (c) 
For partially locked states with K < Ki, there is still an eigenvalue at 0. The rest of the 
spectrum is continuous, as before, but now it includes the whole imaginary axis. Hence 
the partially locked states are linearly neutrally stable. 
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